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1^ I We derive the electron temperature gradient in the Galactic disk using a sam- 

ple of Hll regions that spans Galactocentric distances 0-17 kpc. The electron 
temperature was calculated using high precision radio recombination line and 
^ , continuum observations for more than 100 Hll regions. Nebular Galactocen- 

^ I trie distances were calculated in a consistent manner using the radial velocities 
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measured by our radio recombination line survey. The large number of nebulae 



0> I widely distributed over the Galactic disk together with the uniformity of our 

\Q [ data provide a secure estimate of the present electron temperature gradient in 

O ■ the Milky Way. Because metals are the main coolants in the photoionized gas, 



the electron temperature along the Galactic disk should be directly related to 
the distribution of heavy elements in the Milky Way. Our best estimate of the 



O 

H ' electron temperature gradient is derived from a sample of 76 sources for which 



O ■ 

^ I we have the highest quality data. The present gradient in electron temperature 

has a minimum at the Galactic Center and rises at a rate of 287 ± 46 Kkpc^^. 
^ I There are no significant variations in the value of the gradient as a function of 

' Galactocentric radius or azimuth. The scatter we find in the H ll region electron 

temperatures at a given Galactocentric radius is not due to observational error, 
but rather to intrinsic fluctuations in these temperatures which are almost cer- 
tainly due to fluctuations in the nebular heavy element abundances. Comparing 



^Observatorio Nacional, Rua General Jose Cristino 77, 20921-400, Rio de Janeiro, RJ, Brazil 
(quireza@on.br) . 

^Instituto de Astronomia, Geofi'sica e Ciencias Atmosfericas (lAG), Universidade de Sao Paulo, Rua do 
Matao 1226, 05508-900, Sao Paulo, SP, Brazil. 

^Astronomy Department, University of Virginia, P.O.Box 3818, Charlottesville VA 22903-0818, USA. 

^Institute for Astrophysical Research, Department of Astronomy, Boston University, 725 Commonwealth 
Avenue, Boston MA 02215, USA. 

^National Radio Astronomy Observatory, P.O. Box 2, Green Bank WV 24944, USA. 



- 2 - 



the H II region gradient with the much steeper gradient found for planetary nebu- 
lae suggests that the electron temperature gradient evolves with time, becoming 
flatter as a consequence of the chemical evolution of the Milky Way's disk. 

Subject headings: H II regions — ISM: abundances, clouds, evolution, and lines, 
structure — nucleosynthesis, abundances — radio lines: ISM 

1. INTRODUCTION 

Churchwell & Walmsley (1975) pioneered Hll region radio recombination line (RRL) 
studies of the relationship between nebular electron temperatures, Te, and Galactocentric 
distance, i?gai, (Churchwell et al. 1978; Wink et al. 1983; Shaver et al. 1983; and others). 
Because RRLs are not obscured by interstellar dust, relatively faint Hll regions at extremely 
large distances from the Sun could be detected. They found that there was a Galactic 
temperature gradient wherein Tg is low in the Galactic Center and increases with i?gai. Such 
a gradient was first observed in other nearby spiral galaxies by Searle (1971), Rubin et al. 
(1972) and Smith (1975). This Milky Way electron temperature gradient was confirmed by 
radio continuum emission (Omar et al. 2002) and by [O ill] forbidden line optical observations 
(Peimbert et al. 1978; Deharveng et al. 2000). Because heavy elements cool photoionized 
gas, H II region electron temperatures arc directly related to the heavy element abundance: 
low Te corresponds to higher heavy element abundances because of the greater cooling rate 
and vice versa. Consequently, Tg gradients should be inversely related to metal abundance 
gradients (but see §6). Because there has been more stellar processing in the inner Galaxy 
one expects that on average the metallicity decreases as a function of i?gai. 

The existence of a Galactic gradient in Hll region electron temperature is now firmly 
established. Nevertheless, there are still uncertainties in the magnitude of the gradient and 
the possible existence of real variations, both in it!gai and Galactic azimuth, of nebular Tg. 
Extant RRL studies yield Tg gradients that, roughly, vary from 250-440 Kkpc~^. Discrep- 
ancies in the results obtained from different studies may be attributed to several factors, 
including the source sample, the Te derivation, etc. Since the exact value of the Tg gradi- 
ent provides an important constraint on models for Galactic chemical evolution, it must be 
determined as accurately as possible. 

Here we derive H ll region electron temperatures based on radio recombination line and 
continuum data for a large sample of nebulae widely distributed across the Galactic disk. 
The RRL data are of unprecedented sensitivity compared with previous studies. We exam- 
ine anew the Galactic temperature gradient and assess the magnitude and origin of the Tg 
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dispersion at a given i^gai- In § 2 we describe our H ll region sample. In § 3 we derive the neb- 
ular electron temperatures and discuss non-LTE effects. Nebular Galactocentric distances, 
i?gai, and heliocentric distances, dsun, are derived in §4. In §5 our new determination of 
the electron temperature gradient in the Galactic disk is made. Here we also investigate a 
possible spatial variation of this gradient. We discuss the astrophysical implications of our 
efforts in 5 6. 



2. OBSERVATIONAL SAMPLE 

Our Hll region sample data are described by Quireza et al. (2006a) and references 
therein. The radio recombination line and continuum data result from two different exper- 
iments. Neither was targeted to investigate radial gradients in the Galactic disk. The first 
one (hereafter '^He survey) is a study of the abundance of '^He in the Milky Way interstellar 
medium (Rood et al. 1984; Bania et al. 1987, 1997; Balser et al. 1994) using the hyper- 
fine transition of ^He''' at 8.665 GHz. The observing techniques required that a number 
of radio recombination lines be measured simultaneously with the ^He"*" transition. These 
RRLs were used both to monitor the system performance and also to measure the spectral 
baseline frequency structure. As a consequence of the large integration times (over 100 hours 
in many cases) accumulated during this experiment, wc obtained extremely high sensitivity 
measurements of the He 91 a and H91a (An = 1) RRL transitons for a significant number 
of Galactic Hll regions. The second experiment (hereafter Cll survey) is a study of Cll 
recombination lines in photo-dissociation (PDR) regions surrounding the Hll regions. For 
this survey we simultaneously observed the H, He, and C 91 et and 92 a RRL transitions. 
Since recombination lines of the same order with similar principal quantum numbers, such 
as H91q; and H92q; or He9lQ; and He 92 a, should have the same intensity, we averaged 
the 91 a and 92 a spectra to attain higher sensitivity (see Quireza et al. 2006a for details). 

All observations were made near 8.6 GHz (3.5 cm) with the National Radio Astronomy 
Observatory^ (NRAO) 140 Foot (43 m) telescope in Green Bank, WV, which has a half power 
beam width (HPBW) of 3(20 at this frequency. Our sample has 106 sources: 47 nebulae 
from the ^He survey and 66 from the C ll survey. (There are 7 objects in common.) 

The quantities needed to derive nebular electron temperatures are the H91 o; line peak 
intensities, T^, the line full width at half-maximum, Av, and the continuum intensity, Tq- 
These parameters and their errors may be found in Quireza et al. (2006a). There we also 
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assess the quality of our data: we have defined quahty factors, QFs, for both spectral fine 
and continuum data. Because we cannot quantify the systematic errors, these QFs provide a 
qualitative measure of them. By choosing subsets of data with different QF's we can assess 
the importance of systematic errors on our conclusions reported here. Quality factor A 
sources refer to our best data: the spectra are almost noiseless and systematic errors appear 
to be neglibible. The QF decreases from A to E. We judge QFs D and E to be of too low 
confidence to be included in our analysis here. 



3. DETERMINATION OF THE PHYSICAL PARAMETERS 

3.1. LTE Electron Temperature 

The electron temperature of each nebula is derived from the observed radio recombi- 
nation line-to-continuum ratio, Tl/Tc- We model our sources as homogeneous, isothermal 
spheres. This approximation allows us to calculate the electron temperature without having 
to know the distance to the source. If the distance is known many other nebular physical 
properties, such as electron density and emission measure, can be derived. 

Most of the emission observed from H ll regions is continuum radiation produced by free- 
free thermal Bremsstrahlung in the plasma. At high frequencies the nebular gas is optically 
thin and the ratio between the brightness temperature of a recombination line and that of 

the free-free emission continuum depends of the radio frequency and the gas temperature, 
but is independent of the electron density, rie- Thus the observed Tl/Tc may be used to 
estimate the electron temperature of the Hll region (Goldberg 1968; Rohlfs & Wilson 2000). 
Assuming local thermodynamic equilibrium (LTE) and negligible pressure broadening of the 
lines by electron impacts, the electron temperature is given by: 

0.87 



UJ" '^GHzJ I ri(H+) j V kms-' J I «(H+) 
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where we distinguish the LTE temperature, T*, from the electron temperature corrected for 
non-LTE and high-density effects, Tg. The line frequency, vl = 8.584823 GHz, corresponds 
to the rest frequency of the H 91 a recombination line. Here Tq (K) is the continuum antenna 
temperature; T^ (K) and Av ( km s~^) are the H 91 o; recombination line antenna temperature 
and FWHM hnewidth. The n(^He+)/n(H+) ionic abundance ratio was calculated using the 
areas of Gaussian fits to the H and He recombination lines (Peimbert et al. 1992): 

n(^He+) _ rL('He+) Av{^}ie+) 
n(H+) ~ Tl(R+) Av(H+) ' 
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For a small number of objects we do not have good measurements of the ^He transition. In 
these cases we used a constant value of 0.07 ± 0.02 for the n(^He'^)/n(H+) abundance ratio. 
This ratio is the average of our 80 best QF {A, B, and C) sources and is typical for Hll 
regions (Churchwell et al. 1974; Shaver et al. 1983). 

The nebular LTE electron temperatures and their errors, ar^, together with the Galacto- 
centric distances, -Rgai, heliocentric distances, dgun, and the n(^He"'")/n(H+) ionic abundance 
ratios, are listed in Tabic 1. Also given are the name of the source and some of its physical 
properties including the spherical angular size, ©diami linear diameter, D, flux density, S^, 
continuum brightness temperature, Tq, and electron density, ric, (see the discussion below). 
Each nebula's survey membership (either ^He or C ll) is indicated. 

The electron temperature errors were derived by propagating the Gaussian fltting errors 
for the line and continuum measurements. These (Xt^ errors vary from 0.3-17.8% (2.2% on 
average). For our best QF data, errors vary from 0.4-11.0% (1.3% on average). These 
are lower limits to the temperature errors. Due to baseline problems and complex nebular 
structures, uncertainties in continuum measurements are certainly larger than the one we 
estimate, reaching 10% or even 20% in the worst cases. Since we have no way to quantify 
systematic uncertainties, we use the quality factors to estimate the effect of the systematic 
errors on the electron temperatures. The QFs for continuum and spectral line parameters 
are listed, in this order, in Table 1. 

The continuum observing mode can also affect the uncertainty in the measurements. 
We made continuum measurements using both the switched power (SP) and total power 
(TP) techniques (see Quireza et al. 2006a). There we point out that as a source was tracked 
across the sky we interleaved the recombination line and SP continuum measurements such 
that both the continuum and line data span the same hour angle ranges and experience the 
same weather conditions. Thus the line and continuum data would need identical corrections 
for telescope gain and atmospheric opacity. When calculating source properties that depend 
on the line-to-continuum ratio, such as the electron temperature, Quireza et al. suggest 
using SP continuum measurements. We therefore normally use SP continuum observations 
to calculate T* because many telescope effects are canceled in the hne-to-continuum ratio. 
We did not have SP continuum measiucmcnts for a small number of our sources and had to 
use TP continuum data to calculate T*. We estimate that the use of TP continuum data 
should give at most a ~ 20% uncertainty in T*. The observing mode used is listed in Table 1 
for each nebula. 

Accurate calculations of Tg should incorporate non-LTE effects including departures 
from LTE, stimulated emission, and pressure broadening from electron impacts. These 
corrections for non-LTE effects are sensitive to the local density and thus the Hll region 
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geometry. Nevertheless, it has been shown that under many conditions LTE is a good 
approximation and the LTE electron temperature T* is close to T^.. 

For conditions close to LTE and where pressure broadening is not significant, the (3 /a 
RRL intensity ratio should be about 0.28, and the (3 /a RRL fine width ratio should be close to 
1 (Shaver & Wilson 1979). We determined these ratios for the ^He survey where the H 114/3 
line has been observed. For our best QF data (28 objects), the average intensity ratio is 
0.26 ± 0.03 and the average line width ratio is 1.02 ± 0.06. Moreover, Shaver (1980) defined 
a RRL observing frequency, z^lte, such that T* = Te. This optimal radio recombination 
line frequency is a function of emission measure: z^lte ~ 0.081£'M'^'^^. It is essentially 
independent of density, temperature, or structure within the nebula. Using our best QF 
data (88 objects) we obtain an average frequency z^lte = 6.9 ± 3.6 GHz. This is close to 
our observing frequency of 8.584 GHz. Furthermore, detailed density structure models for 
a subset of our ^He survey nebulae were made using high resolution VLA radio continuum 
images and high-order RRLs that were sensitive to local electron densities (Balser et al. 
1999). The non-LTE electron temperatures calculated for these models are very close to 
the LTE electron temperature determined using the H91q; lines. We therefore conclude 
that our LTE electron temperature must be very close to the real average nebular electron 
temperature. 



We derived the Hll region angular sizes by assuming homogeneous, spherical nebulae. 
Each nebula has an observed full width at half-maximum size of Q^rc = [Q{RA) Q^DEC)]^^"^ 
(geometric mean, see Quireza et al. 2006a). Assuming that the source has a Gaussian 
brightness distribution and that the telescope beam has a Gaussian pattern with half power 
beam width, HPBW, then ©grc is: 



where Qg is the nebular Gaussian angular size (Mczgcr & Henderson 1967). We use the 
technique developed by Panagia & Walmsley (1978) to derive 0diam, the angular size of 
a homogeneous, spherical nebula from Qq- We use TP continuum measurements if at all 
possible. The nebular ©diam so derived is listed in Table 1 together with the finear diameter, 
D. Objects whose angular sizes are close to the telescope HPBW are not well resolved and 
consequently have less precise estimates for their sizes. This will give a larger error for any 
physical property whose derivation depends on the angular size. One important property 
which requires precise angular size measurements is the electron density. 



3.2. Nebular Angular Size 




(3) 
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3.3. R.M.S. Electron Density 

If a homogeneous, spherical nebula is also optically thin, then the electron density can 
be calculated from the peak continuum brightness temperature of the source, Tq (Balser 
1995): 



Tie {cm ^) 



Tc^(K)z/(GHz)2t;(K)1/2 



1/2 



(4) 



8.77 X 10 ^ ©diam (arcmin) dsun ( kpc ] 

where z/ is the frequency (8.66565 GHz), T* is the electron temperature, X = 4.954 x 
10-2 T* '^''^jv, Gdiam is the spherical size, and dsun is the source's distance from the Sun (see 
§6). Here we assume no doubly ionized helium gas within the H ll region. The flux density, 
S^, is: 

{Jy) = 1.223 {K) (^^\ = 2.647 (^] (^] , (5) 



^ J x'nhj \ >^ 

where temperatures are in K, angular sizes are in arcmin, and the wavelength is in cm. The 
beam efficiency of the 140 Foot telescope is r^b — 0.86. 

The nebular flux and electron densities calculated in this way are listed in Table 1. 
Because Tq measures the continuum radiation produced by free-free thermal Bremsstralung 

in the plasma its value depends on the integral of along the line of sight. Thus the electron 
densities derived here are the root mean square (rms) density of the entire Hll region. If 
the nebula is homogeneous, then the rms electron density equals the local density. Detailed 
analysis, however, shows that this is not generally true for Galactic H II regions (Osterbrock 
& Flather 1959). The rms densities are probably somewhat lower than the true densities of 
the emitting regions, since the gas is not evenly distributed throughout the H ll region (Fich 
& Silkey 1991). 



3.4. Compcirison with Previous Studies 

Here we compare our electron temperatures with results from the literature for nebulae 
in our sample. The majority of these previous efforts were also studying the Milky Way 
Te gradient. Most Tg estimates come from RRLs, although measurements in other spectral 
regions are available for a limited number of objects. We summarize the Tg differences in 
Table 2 which lists the number of sources in common, N, the average percentage difference 
between our results and the reference, and the average percentage of the absolute value 
of difference between our results and the reference. The average difference will reveal any 
offset in Tg scales; the average absolute value difference measures the scatter between the 
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studies. The table also lists the type of observation and the method used to calculate Tg 
(line-to-continuum ratio; optical forbidden lines; radio continuum emission; etc). 

It is immediately apparent that our Tg's are systematically high compared with previous 
studies. While one might expect systematic differences between TgS determined by different 
methods, the largest offsets are found between our values and those of Wink et al. (1983) 
(WBW83) and Shaver et al. (1983) (SMNDP83) which are both RRL studies similar to ours 
with many sources in common. In Figure 1 we compare our LTE electron temperatures with 
their Tg 's. Our best QF data {A, B, and C) are plotted as filled symbols. Our values are 
systematically 11% higher than WBW83 and 13% higher than SMNDP83. 

We considered the origin of these differences in some detail. Each and every term 
in Eq. 1 may be the explanation. There may also be issues of calibration between line and 
continuum measurements and between telescopes. All the RRL studies were made at different 
frequencies using different telescopes with correspondingly different beam sizes. Thus each 
survey probes different volumes of each nebula. Certainly there is complex structure inside 
some H ll regions which includes large density and temperature fluctuations. The Tg derived 
from disparate RRL transitions using different telescopes can well be different, especially 
under the assumption of LTE. 

The beam size for the WBW83 sudy was substantially smaller than ours, so on the 
average they were observing higher density gas than we did. Exactly how to interpret this is 
quite complex and would require detailed density and excitation modelling for each nebula. 
The continuum measurement technique is also important. These details can lead to Tq 
measurements that are too low (the continuum data do not extend beyond the Hll region) 
or too high (Galactic non-thermal continuum is included in Tc)- 

In Balser et al. (1999) our continuum observations from the 140 Foot, MPIfR 100 m, 
and the VLA were carefully cross-calibrated, confirming the calibration techniques used for 
the ^He and Cll surveys. Furthermore, our recombination line data are vastly superior to 
these previous studies due to a combination of improved technology (system temperatures of 
~35K compared with ~100K) and longer integration times (tens of hours compared with 
tens of minutes). In sum, we have yet to identify just why WBW83 and SMNDP83 have 
lower Te values. 

We have in fact reason to believe that our electron temperatures are the best values 
derived to date. We have substantially better recombination line data that have unsurpassed 
sensitivity due to a combination of modern receivers and our extremely long integration 
times. The modes of the distributions of source integration times for the carbon and ^He 
surveys are ~ 15 hrs and ~ 50 hrs, respectively. Previous efforts have source integration 
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times of ~ 30 min or less. Finally, these extremely high signal to noise spectra allow us to 
model the spectral baseline frequency structure with unprecidented accuracy which gives us 
great confidence in our determination of the recombination line parameters. 

4. NEBULAR DISTANCES 

We derived the Galactocentric distance, i?gai, for each nebula using the observed re- 
combination line LSR velocity and assuming a Galactic rotation curve. For sources located 
inside the solar orbit we used the Clemens (1985) rotation curve; otherwise the Brand & Blitz 
(1993) rotation curve was used. Both rotation curves assume purely circular rotation laws 
and they place the Sun at a Galactocentric distance of Rq = 8.5 kpc orbiting the Galactic 
Center at an LSR circular velocity of ©o = 220 km s^^. The nebular RRL velocities can 
be found in Table 2 of Quireza et al. (2006a). Here we have thus chosen to use only these 
kinematic determinations of the Rga,i of our sources. Many of the previous studies of the 
Milky Way Tg gradient use a mix of techniques to estabhsh the nebular i?ga,i. 

The nebular heliocentric distance, (isun, is also listed in Table 1. These distances were 
also derived kinematically from the observed nebular LSR velocity. For sources located 
inside the solar orbit, each radial velocity value corresponds to two distances ( the "near" 
and "far" kinematic distances) equally spaced on either side of the tangent point. In most 
cases, we were able to resolve the kinematic distance ambiguity by measuring the 21 cm Hi 
absorption spectrum toward the nebular continuum. Discrimination between the near and 
far distance was done by comparing the maximum velocity of the H I absorption with that 
of the Hll region recombination line. For those objects for which we could not resolve the 
ambiguity via Hi absorption, we used distances available in the literature. Our Hi survey, 
observational technique, method of analysis, and detailed description of our dgun derivations 
are described by Quireza et al. (2006b). 

Figure 2 shows the distribution of our Hll region sample projected onto the Galactic 
plane. Only sources with known rfgun ai'c plotted. The majority of our Hll regions are lo- 
cated in the first and fourth Galactic quadrants and their i?gai's can reach ~19 kpc. Some 
sources are located beyond the Galactic Center with lisun ~ 20 kpc. Different symbols identify 
the ^He (triangles) and C ll survey (circles) nebulae. The ^He sources span a larger range 
of Galactocentric radius and azimuth than do the C ll nebulae. Filled symbols denote our 
highest QF sources (QFs A, 5, & C); open symbols flag our poorer quality data. Unfortu- 
nately, many of the low QF objects are astrophysically significant; they are located in areas 
of the disk not well covered by nebulae with more accurate data. 
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There is good reason for using only kinematic distances: homogeneity of approach. A 
detailed comparison between optical and kinematic distances is beyond the scope of this 
paper. Both approaches have their limitations and perhaps the resolution will be a direct 
trignometic measurement using VLBI techniques which is being pioneered by Mark Reid and 
his team (Xu et al. 2006; Hachisuka et al. 2006). 

It is common to cite optical spectrophotomctric distance uncertainties to be of order 
15%. But as one delves into the literature in detail one finds a greater dispersion. For 
any given analysis the 15% holds, but slight differences in spectral classification, adopted 
luminosity scale which is model atmosphere dependent, and uncertainties in the extinction 
correction conspire to give a larger uncertainty when different research groups measure the 
distances to the same object. These systematic errors are the optical counterpart of the 
streaming motion uncertainty inherent in kinematic distances. So the consistency of optical 
distances may not be all that much better than the dispersion between optical and radio 
distances. Most of our H II regions at not optically visible. A quick web search for the 7 most 
distant (in Rgax) nebulae shows that the mean absolute value of the radio-optical discrepancy 
for these sources is 11% ± 6%. 

Non-circular streaming motions will certainly affect the kinematically derived distances. 

If we knew what the Galactic scale streaming motions arc wc would certainly correct for them. 
The Galactic Bar, the spiral density wave, accretion events such as the Sgr dwarf elliptical 
galaxy, and close encounters such as that which occurcd for the LMC will all generate large- 
scale, asymmetric streaming motions in the Milky Way disk. At present there is no way to 
calculate these effects accurately and thus generate the true velocity field of the Galaxy. 

5. ELECTRON TEMPERATURE GRADIENT 

5.1. Gradient Magnitude 

Our H II region sample contains some sources observed in both the ^He and C ll surveys 
as well as nebulae that have RRL spectra taken toward multiple positions. For sources 
common to both surveys^, the T* and Rga,i values derived for the ^He and Cll surveys were 
averaged (simple mean) . These nebulae do not appear twice in the analysis of the electron 
temperature gradient. We only averaged data of good QFs (Cor better) and which also were 
in good agreement. (The origin of some differences between the two surveys is discussed by 



2G16.936+0.75 (M16), G23.421-0.21, G25.382-0.17 (3C385), G43.169-F0.0 (W49), G49.384-0.29, 
G49.486-0.38 (W51) and G79.293+1.29 
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Quireza et al. [2006a].) These T* values are consistent with each other within the errors; 
differences are not larger than ~ 10%. 

In the ^He survey five morphologically complex H ll regions had spectra taken toward two 
or three different positions. These nebulae are: G16. 936+0. 75 (the brightest and northern 
most components: M16 and M16N), M17 (northern and southern most components: M17N 
and M17S), NGC6334 (NGC6334A and NGC6334D), Rosette (Rosette A and Rosette B) 
and S209 (brightest, northern and southern most components: S209, S209N and S209S). 
Because these nebulae are extended objects with angular diameters larger than the 3^20 
HPBW beam size, each of the observed components is included in our analysis of the T* 
gradient. Since these components are separated by more than a beamwidth any T* differences 
within an object are real temperature fluctuations. 

Figure 3 shows our nebular LTE electron temperatures, T*, plotted as a function of the 
Galactocentric distance, i?gai. We include only the 78 sources with our best data (QFs = A, 
B, h C — Sample B described below) for both line and continuum. Shown are least squares 
linear fits to the gradient, T* {K) =01 + 02 -Rgai (kpc), for the entire sample and, separately, 
for nebulae located inside and outside the solar orbit. The gradient is flatter for H ll regions 
in the inner Milky Way. 

Seven nebulae in Figure 3 are flagged because they do not follow the general T* gradient. 
Two nebulae, G49.582-0.38 (1,851 K, 6.5 kpc) and G5.956-1.265 (3,416 K, 7.8 kpc) have 
temperatures much lower than the 7, 585 ± 1, 262 K average T* for the 6-8 kpc interval of 
Rgai- Moreover, the G49. 582— 0.38 temperature is ~ 3,600 K lower than the Wink et al. 
(1983) value. For these sources we had to use TP continuum measurements to derive T*. 
This may have compromised the accuracy of our result. G5. 899— 0.427 (11,128 K, 6.1 kpc) 
also lies far from the general trend of the sample. Here, however, we do not have any reason 
to suspect that its T* derivation is less reliable than that of other sources in the 6-8 kpc 
zone. 

The four remaining anomalous nebulae are located at the extremes of the i?gai distribu- 
tion of our H II region sample. These nebulae thus have a large influence on the temperature 
gradient fits. The two Galactic Center sources, Gl. 13—0.1 (7,135 K, 0.1 kpc) and Sgr B2 
(8,169 K, 0.4 kpc), may well share the anomalous chemical abundances in this region (see 
§6.2). The two outer Galaxy sources, S209 (10,506 K, 16.9 kpc) and S209 N (12,565 K, 16.7 
kpc), are part of the same Hll region. Here again we can identify no compelling reason to 
exclude these nebulae. 

We analyzed 6 subsets of our H ll region sample in order to assess the effects of sample 
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selection on the temperature gradient.^ Sample A includes all 109 of our sources. Sample B 
includes only high QF sources (the 78 sources shown in Figure 3 which have QFs C or better). 
Sample C removes TP continuum sources from sample B for a total of 64 nebulae. Sample D 
removes the 7 sources flagged in Figure 3 from Sample B for a total of 71 nebulae. Sample E 
removes G49. 582— 0.38 and G5. 956— 1.265 from sample B for a total of 76 nebulae. Sample 
F removes 3 sources from sample E (the two Galactic Center sources and G5. 899— 0.427) for 
a total of 73 nebulae. 

We fit least squares linear temperature gradients to each sample. The fit results are 
summarized in Table 3 which lists the coefficients Oi and 02 and their the standard deviations 
(j{ai) and £1(02), together with the correlation coefficient, r, the of the fit, the number, A^, 
of sources in the sample, and the i?gai range of the sample, A i?gai. In all cases the gradient 
is close to 300 Kkpc~^. Even our smallest sample of 64 objects is astrophysically significant 
in the sense that it spans a large range of Galactic radius. The fit to the sample E nebulae 
is shown as a solid line in Figure 3. Unless otherwise stated we use the sample E gradient 
fit in all subsequent analysis described herein. This fit is Tg = (5780 ± 350) + (287 ± 46)i?gai 
(r = 0.59, Ar = 76). 

Figure 3 shows T* error bars. Many nebulae have error bars which are smaller than the 
plotted symbols. The error bars shown arc the propagated statistical errors in the measured 
quantities. Various systematic effects almost certainly lead to uncertainties larger than the 
error bars shown. Because systematic errors are inherently unquantifiable and certainly are 
not normally distributed, any attempt to plot "systematic error bars" would be misleading. 

This is especially true for errors in the kinematically determined i?gai. They are entirely 
dominated by systematic effects including Galactic scale streaming motions, differences in 
the choice of rotation curve, etc. Estimates for the magnitude of the error of the i?gai 
determination are typically ~25% (e.g. Kuchar & Bania 1990). This value is the quadrature 
sum of the ~ 15% error due to the uncertainty in the rotation curve and of the ~ 20% error 
due to streaming/non-circular motions. Here we have, however, at least made the i?gai 
determinations in a uniform, systematic way that distinguishes this effort from the majority 
of Tg gradient studies which normally draw their distances from a heterogeneous mixture of 
techniques. 

It is important to understand the nature of the scatter in Figure 3. Deharveng et al. 
(2000) argued that a significant fraction of the scatter in earher investigations of the electron 



^G345.40+1.41, S291 and S132 were excluded from all samples. We have no continuum information for 
G345.40+1.41. The other nebulae have large uncertainties associated with either their line or continuum 
data. 
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temperature gradient may result from observational errors and sample inhomogeneity. They 
cite the case of the bright H ll region S206: the electron temperature estimated by various 
authors ranges from < 8, 000 K (Churchwell et al. 1978; Mezger et al. 1979) to about 13,000 
K (Litchen et al. 1979). In our study sample homogeneity is not a problem and the statistcal 
errors are quite small compared to the scatter in T*. 

Could systematic error be responsible for the scatter? We believe that the main source of 
systematic error arises from determining the baseline level of our continuum measurements. 
The continuum QF's givem at least a qualitative estimate of this. In Figure 4 we plot the 
Figure 3 points without error bars and use different symbols to identify the continuum QFs. 
The two lower outliers are indeed QF = C, but for the bulk of the points the differing QF's 
yield a comparable scatter. From this we conclude that the scatter in the points indicates a 
real spread in T*. That is, the dispersion in T* at any i?ga,i is indicative of intrinsic variations 
in T* between nebulae. 

5.2. Electron Temperature Vetriations in Galactic Radius 

Given the size, uniformity, and precision of our sample, we can investigate whether there 
is more complex behavior than a linear gradient in the variation of T* with i?gai. To reduce 
the scatter to make trends more visible we have "smoothed" the data in two ways as shown 
in Figure 5. The symbols with error bars are a 10 point running mean of T* and it!gai plotted 

at intervals of 4 points along Galactocentric distance. The horizontal line segments show the 
mean T* in each 1 kpc interval, ranging from i?gai = 0-1 kpc to i?gai = 16-17 kpc. These are 
offset by +2,000 K for clarity. Both techniques smooth the data in slightly different ways; 
both suggest a slightly smaller gradient in the inner area of the Galactic disk. 

We also made a least squares second order polynomial fit to the electron temperature 
gradient for the sample E sources. This fit was very similar to the best linear fit, and perhaps 
even misleading because of the effect of the radial outliers. We therefore divided sample E 
sources into nebulae inside and outside the solar orbit and then fit two linear segments to 
these data subsets. These fits are shown in Figure 3. The slope interior to i?gai is less than 
that in the outer Galaxy: 153 ± 85 Kkpc~^ compared to 404 ± 130 Kkpc~^. This result 
is strongly influenced by the outliers. If we exclude Gl. 13-0.1, G5.899-0.427, Sgr B2 and 
S209 from the hnear fits, we have a gradient of 268 ±96 Kkpc"^ for i?gai < Rq and 342 ±239 
Kkpc~^ for i?gai > -Ro- We conclude that our data do not justify anything more elaborate 
than a single linear fit. The data do hint at more complex behavior, but a larger sample 
of nebulae at both very small and very large Rgsi are needed to explore it. This is to be 
expected since our nebular sample was not chosen to study the disk electron temperature 
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gradient. Indeed having more nebulae at larger Galactocentric distances would significantly 
improve the determination of the electron temperature gradient. We in fact intend to make 
observations of more nebulae at larger Rggx in the future because of this. 

5.3. Electron Temperature Variations in Galactic Azimuth 

In § 5.1 we conclude that the scatter in T* at a given i?gai is not due to observational 
error. Figure 6 shows a histogram of the percentage deviation of the nebular T* from the best 

Te gradient model fit to sample E. The deviations are Gaussian distributed with a ATg/Te 
dispersion of about 14%. This implies that Galactic H II regions have intrinsic T* fluctuations 
of ~ 1,100 K at any -Rgai- Shaver et al. (1983) also argued that most of their scatter in 
is intrinsic. Their Figure 17 shows that a realistic range in the effective temperature of the 
exciting stars (30,000-45,000 K) or of the electron density can account for a spread of Tg as 
large as 2,000 K; this is confirmed by photoionization models (Rubin 1985; see §6.1). 

Maps of radio continuum emission (Altenhoff et al. 1978; Reich et al. 1990; Fiirst et al. 
1990; Haynes et al. 1978), and our own continuum observations show that many H II regions 
are found spatially close to each other. The gas within a given complex of nebulae could share 
the same nucleosynthetic history perhaps including self-enrichment. The nucleosynthetic 
history might, however, vary significantly from complex to complex. We searched our Hll 
region sample for spatially clumped clusters of nebulae with similar T* values. We found no 
obvious signature of patchy nucleosynthesis. 

Another possibility is that the assumption of axial symmetry in the stellar production of 
heavy elements is invalid. If the radial gradient were a function of Galactocentric azimuth, 
then the scatter in a T* vs i?gai Figure 3-type plot would result from lumping together 
nebulae from different azimuths into the same i?gai bin. We have therefore searched for 
azimuthal differences in our nebular sample. 

Figure 7 shows the nebular electron temperature plotted as a function of Galactocentric 
radius for four distinct ranges of Galactocentric azimuth. A difficulty in the analysis of 
the azimuthal variation of the temperature gradient is that our sample is not uniformly 
distributed in the disk. Some azimuth intervals have large concentrations of Hll regions 
(e.g. 350°-20°), while others have no objects at all (see Figure 2). We divided our Hll 
region sample into four azimuth ranges chosen such that each contains a comparable number 
of nebulae: 300°-30°, 30°-90°, 90°-150°, and 150°-215° (G49.582-0.38 and G5.956-1.265 
are not included in this analysis). 

Most of the intervals are relatively well populated between Rg^x 2-10 kpc. Because of 
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this we fit the T* gradient over a shorter i?gai interval, roughly from 3-9 kpc. Properties 
of the gradient fits are given in Table 4 for each azimuth interval. (Table 4 gives the same 
fit information as Table 3.) There is a variation in the T* gradient derived for the different 
azimuth ranges. It is impossible, however, for us to draw any significant astrophysical in- 
ferences from this result because the azimuth ranges do span different angular zones, cover 
different i?gai ranges, and contain different numbers of nebulae. 

In sum, we find no definitive evidence for dumpiness or azimuthal variations in the 
distribution of nebular electron temperatures. A much larger H II region sample that is more 
uniformly distributed in the Milky Way disk is needed. 

5.4. Compcirison With Other Milky Way Electron Gradient Studies 

We searched the literature for previous determinations of the electron temperature gra- 
dient in the Galactic disk. Table 5 summarizes the results of these efforts. Besides the 
reference for and gradient found by each study. Table 5 lists the Galactocentric distance 
interval together with the assumed radius of the solar orbit, Ai?gai(kpc) (Rq [kpc]), the 
number, N , of objects, as well as the type of observation and analysis method used to de- 
rive the gradient. Most studies of gradients in electron temperature in the Milky Way are 
based on observations of radio recombination fines. These RRL gradients vary from 250-440 
Kkpc~^. Wink et al. (1983) suggest a gradient of 270 Kkpc~^ based on observations of the 
76 a line for 84 H ll regions. 

Shaver et al. (1983) derived electron temperature and abundance gradients which have 
often been used as the basis of many models for the chemical evolution of our Galaxy (Tosi 
1988; Giovagnoli & Tosi 1995; Thon & Meusinger 1998; among others). Because their 
gradient, 433 ± 40 Kkpc~^, is one of the steepest values ever determined, here we try to 
understand why this is so. 

Different adopted distances are one possible source of disageeement. The Shaver et al. 
distances were derived from the measured radial velocities using the Schmidt (1965) rotation 
curve with Rq = 10 kpc. According to the authors their Rg^i distances should be accurate 
to within 1-2 kpc in most cases. There are a few objects in common with our sample 
which have i?gai differences much larger than that expected from the change in Galactic 
rotation parameters and rotation curves. For example, G1.13— 0.1 is placed at i?gai ~ 6.0 
kpc in Shaver et al. whereas we derive a distance of i?gai ~ 0.1 kpc. Rudolph et al. (1997) 
give i?gai ~ 0.91 kpc for this Hll region (Lis 1991, Simpson & Rubin 1990). Distances for 
GO. 6— 0.6 are also quite different. One way to compare our results to Shaver et al.'s is to 
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consider only the group of 22 H ll regions in common. If we adopt their distances for these 
22 sources our electron temperature gradient steepens from 258 ± 50 Kkpc"-*^ to 363 ± 37 
Kkpc~^ which is still far short of the 433 ± 40 Kkpc~^ Shaver et al. result. 

The Shaver et al. source sample consists of 67 distinct Galactic H ll regions located in 
the range Rgax — 3.5-13.7 kpc. The RRL data for 44 Hll regions were their own 5 GHz 
(H109q;) observations. This was supplemented by 14.7 GHz (H76q;) observations of 23 
southern Hll regions by McGee & Newton (1981). All these observations were made with 
the Parkes 64 m radio telescope. 

We compare our T* results with those of Shaver et al. in Figure 8 where temperatures 
derived from the H 76q; and H 109a; lines are plotted with different symbols (triangles and 
circles, respectively). To avoid the distance issue, we used their measured velocities to derive 
new kinematic values for i?gai in exactly the same way that we derived our own distances. 
We show only those points for which we could recompute the Shaver et al. Rgai. Our new 
distances (as well as most of the original Shaver, et al. distances) were derived using Galactic 
rotation curves derived from northern hemisphere data. North-South symmetries in Galactic 
rotation have long been known (e.g.. Burton 1988), so it might not be appropriate to use 
our adopted rotation curve for distances to fourth quadrant Hll regions. Fourth quadrant 
sources make up most of the Shaver et al. sample for it!gai > 6 kpc. 

Figure 8 shows that the electron temperatures obtained using the H 76q; data are com- 
pletely consistent with our results. The discordance between our work and that of Shaver et 
al. (1983) arises mostly from the low temperature points calculated with the H 109a lines 

located between Rgsi = 4-6 kpc. Deharveng et al. (2000) also noted that the Shaver et 
al. (1983) electron temperatures in the /2gai = 3-7 kpc zone are generally lower than those 
derived by other RRL studies (Mezger et al. 1979; Wink et al. 1983; Caswell & Haynes 
1987). 

There are at least three factors which could steepen the Shaver et al. (1983) gradient: 
a systematic difference between the Te derived from the H 76a; and H 109a lines; a change 
in distance scales at roughly Rgg^ > 6 kpc; or a different abundance gradient for Fourth 
Quadrant H ll regions."^ 

We also compare our results with those by Afflerbach et al. (1996), who measured 
electron temperatures in ultracompact Hll regions, using the H42a, H66a, H76a and 
H93a RRLs. They found a Galactoccntric gradient [K] = (320 ±64) R^^i + (5537 ±387) 
which is in good agreement with our value despite the fact that ultracompact Hll regions 



''We find an indication of just this for our Fourtfi Quadrant nebuiar sample sliown in Figure 7. 
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are much denser nebulae than our sources. 

Optical studies of Galactic scale H ll region electron temperature gradients are difficult 
because dust in the Galactic plane causes high extinction in the visible. Optical observations 
by Peimbert et al. (1978) gave a large gradient of ~ 1100 K kpc~^. Their result was probably 
influenced by the short i?gai interval spanned by their sample. A more recent optical estimate 
by Deharveng et al. (2000) yields a gradient consistent with that gotten by RRL methods. 
Finally, radio continuum observations at low radio frequency by Omar et al. (2002) are 
also consistent with the gradient estimated by Deharveng et al. (2000) for the i?gai interval 
spanning 10-18 kpc (albeit for a much smaller sample). 

Maciel & Faundez-Abans (1985) investigated the radial electron temperature gradient 
for a large sample of Peimbert (1978) Type II planetary nebulae (PNe) using electron tem- 
peratures derived from forbidden lines of [Olll]. Because Type II PNe have approximately 
circular orbits and thus do not appreciably change their Galactocentric distances during 
their lifetimes, they are well suited to abundance gradient studies. The observed scatter in 
Tg is larger than for H ll regions, probably due to: the large range of effective temperatures 
of the central stars, winds, optical depth effects, etc. Nonetheless, Maciel & Faundez-Abans 
(1985) found a correlation between the electron temperature and Galactocentric distance for 
type II PNe with a gradient of the order of 600 Kkpc~^ and an uncertainty of about 20%. 

Compared to these much older PNe, Hll regions are of zero age and thus sample the 
physical state of the current interstellar medium. Thus the flattening of the vs i?gai 
gradient seen when one compares the PNe and Hll regions may be caused by Galactic 
scale temporal chemical evolution. This time flattening of the electron temperature gradient 
should be accompanied by a corresponding time flattening of the abundance gradient. In fact, 
such a flattening based on the comparison of data from several types of objects (planetary 
nebulae, Hll regions, open clusters, cepheids and young stars) is proposed by Maciel et al. 
(2003; 2005; 2006). It is also predicted by some recent inside-out formation scenario chemical 
evolution models (Hou et al. 2000; Ahbes et al. 2001). 

6. DISCUSSION 

6.1. The Electron Temperature of Galactic HII Regions 

The electron temperature of an H ll region in thermal equilibrium is set by the balance 
of competing heating and cooling mechanisms. It is therefore somewhat surprising that there 
is a Galactic gradient. There are at least four physical properties that could effect Te- (1) 
the effective temperature of the ionizing star, T^Si which sets the hardness of the radiation 
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field exciting (and heating) the nebula; (2) the electron density — coUisional de-excitation 
in the high rig Hll region will inhibit cooling and increase Tf,; (3) dust grains which effect 
the heating and cooling in complex ways; and (4) heavy element abundance which increases 
cooling and decreases Tg (Garay & Rodriguez 1983). 

Rubin (1985) explored how metaUicity, gas density, and the stellar effective temperature 
effect the average electron temperature in model H ll regions. These models predict changes 
in Te of 7,000 K for a factor of 10 change in metal abundance, 2,900 K for a change in 
density from 100 to 105cm-'\ and 1,300 K for a change in T^g from 33,000 to 45,000 K (BO 
to 05 spectral type). Dust grains arc known to play a significant role in the heating and 
cooling of Hll regions (e.g., Mathis 1986; Baldwin et al. 1991; Shields & Kennicutt 1995). 
Photoelectric heating occurs as electrons are ejected from dust grains while the gas is cooled 
by the collisions of fast particles with grains. The electron temperature will decrease with 
distance from the star as the ionizing radiation field is attenuated by dust grains. But the 
electron temperature will also increase as coolants are depleted onto dust grains. Taking 
these competing factors into account Oliveira & Macicl (1986) conclude that dust grains do 
not significantly contribute to the observed Tg gradient with a maximum variation of 500 K. 
Therefore metallicity is the most sensitive factor that sets the nebular T^. value and so Hll 
region metal abundance variations arc the best interpretation for the observed gradient. 

The electron temperature need not be a constant, however, and the observational meth- 
ods used to determine arc sensitive to different regions of the nebula. The two main 
methods used to determine electron temperatures in Hll regions are (1) recombination line- 
to-continuum ratios, such as hydrogen RRL and continuum emission; and (2) forbidden line 
ratios, such as the ratio [O ill] A 4363 / (A 4959 -|- A 5007). The recombination line method 
is weighted towards the lower temperature regions with a weak dependence on Tg, while 
the forbidden line method is weighted towards the higher temperature regions with a strong 
dependence on Tg (e.g., Peimbert 1967). Therefore if temperature structure exists in Hll 
regions these methods can produce different electron temperatures. Observations of [O ill] 
and oxygen recombination lines in Orion show such temperature fluctuations (Esteban et al. 
1998). 

Comparison of hydrogen RRL and [Olll] electron temperatures is more complicated 
since the emission lines arise from different species. The intensities of the RRLs of H relative 
to the underlying continuum, for example, should give an estimate of the electron tem- 
perature in the whole ionized H region (H"*" zone), while observations of the [Olll] optical 
emission should give an estimate of the electron temperature in the O^"*" zone. Photoion- 
ization models indicate that the electron temperature should be higher in the outer regions 
since the ionizing radiation becomes harder with distance from the exciting star and the 
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very efficient coolants like 0+''" are located close to the exciting star (Stasinska 1980; Gar- 
nett 1992). This suggests that Te(0++) < Te(H+) < re(0+) (Stasinska 1990; Stasinska & 
Shaerer 1997; Deharveng et al. 2000). The temperature measured via the RRL, however, is 
not strictly Te(H+) and that measured with the optical [Olll] forbidden lines is not strictly 
Te{0~^~^). Because the recombination lines are emitted preferentially in low-temperature 
regions, Tg (radio) < Te(H"'"), and because the optical forbidden lines are enhanced in high 
temperature regions, Te([Olll])> Tg{0~^~^). Despite this complexity the electron tempera- 
tures determined by these two methods towards the same Hll region produce values of 
that are similar to within the uncertainties (Deharveng et al. 2000). 

Wink et al. (1983) find no correlation between Tg and either the Lyman continuum 
flux, A^L, used to probe different stellar effective temperatures, or the electron density, fig. 
Theoretically, higher values of A^l and Ue should increase (Rubin 1985). But there are Hll 
regions with high A^l and low Ue (e.g., the Rosette nebula) that will increase the scatter in any 
such analysis. Shaver et al. (1983) found a correlation with by dividing the Hll regions 
into two groups of either high or low values of both TVl and Ug- But the correlation seems 
weak since many of the points at smaller Rg^i where the correlation is best arc nebulae 
that appear to have systematically lower electron temperatures (i.e. the H 109 a survey 
sources). The ultracompact Hll region survey of Afflerbach et al. (1996) probes nebula of 
much higher electron density than other surveys. Because these observations were made 
with an interferometer any diffuse gas within their primary beam will be spatially filtered. 
Their derived electron temperatures are higher by ~ 1,000 K compared with other surveys 
of classical H ll regions. 

We explored the effects of A'l and rig on electron temperature in our sample and find no 
correlation. We used an approach similar to Shaver et al. (1983) wherein we selected sources 
either below or above a threshold value of the excitation and electron density (log(NL) = 
49.5 s~^ and rig = 150 cm~^). RRL and continuum surveys of classical H ll regions with single- 
dish telescopes probe Hll region complexes that typically contain more than one ionizing 
star with different spectral types. Such regions consist of both low and high density gas (e.g., 
Balser et al. 1999). But does this account for the observed dispersion of the Tg gradient as 
has been suggested (e.g., Shaver et al 1983)? 

Another possibility is that the dispersion is caused by variations in metallicity. After 
all, stellar abundances at a given i?gai show a dispersion in abundance. Interpretation of 
the stellar abundances is complicated because samples typically contain stars of varying age 
which may well have been formed at a different i?gai than their current location (Edvardsson 
et al. 1993). Correcting for birth location and age Edvardsson et al. (1993) find a scatter 
of 0.5 dex in the solar neighborhood. Priel et al. (2002) find a scatter of about half that in 
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disk star cluster abundances. A scatter of 0.3 dex in heavy element abundance will produce 
a scatter of 2,000 K in Tg, entirely consistent with our observed scatter. 

For Hll regions, Mehringer et al. (1993) suggest differences in metaUicity as the best 
candidate for the ~ 2, 000 K change in Tg between components 1 and 2 in the Sgr B complex. 
Since our survey spans a large region of the Galaxy, both radially and azimuthally, we were 
able to explore such effects by comparing Tg at a given Rg^i over a large range of azimuth. 
The results were inconclusive. Sensitive RRL observations in critial areas of the Galaxy 
could reveal such variations and place important constraints on Galactic chemical evolution 
models. 



6.2. The Galactic Center 

Our nebular sample F fit to the electron temperature gradient excludes the Galactic 
Center (GC) sources SgrB2 and G1.13— 0.1 and includes the outer Galaxy nebulae S209 and 
S209 N. This sample produces the steepest T* gradient, thus providing an extreme constraint 
on Galactic chemical evolution models. 

Early RRL observations of H ll regions in the GC produced electron temperatures that 
were higher than those expected from an extrapolation of the electron temperature gradient 
in the Galactic disk (Mezger et al. 1979; Downes et al. 1980; Wink et al. 1983). Since 
the chemical evolution of the GC may be different than the disk, it is not surprising that 
electron temperatures in the GC do not follow the gradient in the disk. Higher spatial 
resolution observations reveal a wide range of electron temperatures (~ 4, 000 — 10, 000 K) 
in Galactic Center H ll regions (Roelfsema et al. 1987; Gaume & Claussen 1990; Mehringer 
et al. 1992,1993; De Free et al. 1995, 1996; Lang et al. 1997, 2001). The large range of 
electron temperatures may be due to the complexity of the GC. Real variations in mctallicity, 
electron density, and excitation will effect Tg. Some of the higher electron temperatures 
arc overestimated when the radio continuum is contaminated with non-thermal emission. 
Nevertheless, in some objects very high electron temperatures (Tg ~ 20, 000 K) appear to be 
real where the gas may not be in equilibrium (Mehringer et al. 1995). Since the zero level 
for our radio continuum temperature was determined away from the thermal object but not 
far from the Galactic disk, any smooth non-thermal emission should have been removed. 
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6.3. The Galactic O/H Abundance Gradient 

The electron temperature may be converted to oxygen abundance using a correlation 
between metal abundance and temperature. Because of the broad interest in abundance 
gradients, we provide here one such conversion of our T* gradient to an O/H gradient, 
deferring a more extended discussion to a future paper. Using the relation between O/H and 
Te of Shaver et al. (1983) yields the result shown in Figure 9. The O/H gradient is about 
—0.04 dex kpc~^ in agreement with results from Afflerbach et al. (1996), Deharveng et al. 
(2000), Pilyugin (2003) and Esteban et al. (2005). The exact fit to the 76 sample E nebulae 
in Figure 9 is: 

log (O/H) + 12 = (8.958 ± 0.052) + (-0.043 ± 0.007) Rg^i (kpc) (6) 

with a correlation coefficient of r = —0.59. If our observed T* dispersion of ~ 1,100 K is 
caused by real abundance fluctuations, then this would correspond to an [O/H] abundance 
dispersion of ~0.16 dex. 

7. SUMMARY 

We used extremely sensitive radio recombination line observations of a large sample 
of Galactic Hll regions to study the nebular electron temperature gradient in the Milky 
Way disk. LTE electron temperatures were derived from the observed H91a RRL line-to- 
continuum intensity ratios. Departures from LTE were found to be small given the properties 
of our H II regions and the fact that at 8.6 GHz the effect of pressure broadening is not large. 
Our 109 source sample consists mostly of classical H ll regions of low density and with angular 
sizes larger than the 3^20 telescope beam. 

Our temperature gradient analysis has the following virtues: (1) We analyze a large 
sample of Hll regions spanning the entire Galactic disk from i?gai = 0-17 kpc. (2) We use 
high sensitivity radio recombination line measurements obtained with very long integration 
times for our derivation of the electron temperature. (3) All nebular Galactocentric distances 
are calculated kinematically in a self-consistent way. (4) All of the data were observed with 
the same telescope, the NRAO 140 Foot, identically cahbrated, and analyzed in a uniform, 
self-consistent way. 

Our best estimate of the electron temperature gradient is derived from a sample of 76 
sources for which we have the highest quality data (sample E of Table 3). We conclude 
that the present gradient in electron temperature in the Galactic disk has a minimum at 
the Galactic Center and rises at a rate of 287 ± 46 Kkpc~^. This value is consistent with 
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determinations in the optical by Deharveng et al. (2000) and in the radio by Wink et al. 
(1983), Afflerbach et al. 1996 and Azcarate et al. 1985. Our gradient is about a factor of 
two less steep than the gradient proposed by Shaver et al. (1983). 

We find little if any variation of the electron temperature gradient with Galactocentric 
distance. There is some variation of the temperature gradient calculated for different re- 
gions of Galactocentric azimuth. Unfortunately our nebular sample is not homogeneously 
distributed in the Galactic plane which complicates this analysis so no firm conclusions can 
be drawn. 

The scatter we find in the H ll region electron temperatures at a given Galactocentric 
radius is not due to observational error, but rather to intrinsic fluctuations in these tempera- 
tures that are almost certainly due to fluctuations in the nebular heavy element abundances. 
A comparison of the Hll region gradient with much steeper gradient found for planetary 
nebulae suggests that the electron temperature gradient evolves with time, becoming flatter 
as a consequence of the chemical evolution of the Milky Way's disk. 
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9120 ± 40 


SP 


A A 


Cii 


G93.060+2.81 




12.3 


8.5 






0.068 ± 0.024 


0.0 






10840 ± 270 


TP 


C C 


3 He 


G102.88-0.72 


S132 


11.3 


5.8 


17.0 


29.0 


0.075 ± 0.035 


7.8 


0.3 


17.3 


14800 ± 12780 


TP 


E E 


3 He 


G107.18-0.95 


S142 


10.6 


4.3 


10.7 


13.4 


0.052 ± 0.018 


4.6 


0.4 


29.1 


10400 ± 270 


TP 


D E 


3 He 


GllO.11+0.04 


S156 


11.5 


5.3 


1.8 


2.8 


0.040 ± 0.010 


1.7 


5.4 


232.6 


9070 ± 140 


SP 


A D 


3 He 



Table 1 — Continued 



Source 


Name 


(Kpc) 


Ltsun 

(kpc) 


^diam 
(. ) 


L) 


7i( tie )/n(ti ) 


a 




Tie 

(cm^-^) 




iVLoae 




Survey 


Gill. 53+0. 82 


S158/NGC 7538 


12.4 


6.4 


4.8 


8.8 


0.089 ± 0.004 


18.9 


8.2 


158.2 


8230 ± 40 


SP 


B A 


3 He 


G112.24+0.23 


S162/NGC 7635 


11.1 


4.7 


6.8 


9.2 


0.083 ± 0.015 


3.4 


0.7 


46.1 


8070 + 130 


SP 


B C 


■'He 


G118.15+4.96 


S171 A 


9.0 


0.9 


22.8 


6.1 


0.064 ± 0.012 


24.7 


0.5 


46.4 


9540 ± 1050 


TP 


C C 


•'He 


G133.72+1.21 


W3 


11.7 


4.1 


3.1 


3.7 


0.075 ± 0.004 


40.0 


41.1 


549.0 


8380 ± 40 


SP 


B A 


•'He 


G133.790+1.4 




12.5 


5.0 


4.7 


6.9 


0.081 ± 0.011 


14.1 


6.2 


157.1 


8880 ± 120 


SP 


C A 


3 He 


G150.59-0.95 


S206 


11.5 


3.3 


4.8 


4.5 


0.092 ± 0.008 


5.8 


2.5 


125.0 


9710 ± 90 


SP 


B B 


3 He 


G151. 587-0. 2=^ 


S209 N 


16.7 


8.7 


5.6 


14.1 


0.081 ± 0.019 


7.8 


2.5 


73.7 


12570 ± 360 


SP 


C C 


3 He 


G151. 59-0.23"^ 


S209 


16.9 


8.9 


5.6 


14.4 


0.078 ± 0.007 


7.8 


2.5 


70.7 


10510 ± 90 


SP 


B B 


3 He 


G151. 636-0.5'* 


S209 S 


19.9 


12.0 


5.6 


19.5 


0.079 ± 0.028 


1.0 


0.3 


21.6 


8680 ± 340 


SP 


B E 


•^He 


G155.36+2.61 


S212 


16.7 


8.6 


5.8 


14.4 


0.117 ± 0.023 


2.1 


0.6 


35.5 


10460 ± 240 


SP 


B D 


•'He 


G169. 19-0.90 


S228 


13.8 


5.3 


4.8 


7.4 


0.065 ± 0.028 


1.0 


0.4 


40.3 


9700 ± 740 


SP 


B E 


•'He 


G173.60+2.80 


S235 


10.1 


1.6 


6.6 


3.1 


0.049 ± 0.010 


3.3 


0.7 


81.6 


8940 ±170 


SP 


A D 


3 He 


G189.97+0.40 


S252 


12.1 


3.6 


12.1 


12.8 


0.055 ± 0.017 


5.2 


0.3 


27.6 


9460 ± 230 


TP 


C E 


3 He 


G201.6+1.6 




12.7 


4.4 


10.9 


14.0 


0.078 ± 0.021 


4.1 


0.3 


26.1 


9140 ± 220 


SP 


C D 


3Hc I 


G206. 122-2.3'* 


Rosette B 


10.7 


2.4 


10.0 


6.9 


0.061 ± 0.019 


2.8 


0.3 


32.3 


7480 ± 210 


TP 


D D 


3Hc to 

■^Hc 7 


G206.618-1.8'* 


Rosette A 


10.5 


2.2 


18.0 


11.6 


0.069 ± 0.021 


14.0 


0.4 


32.9 


10770 ± 310 


TP 


D D 


G209.01-19.4 


Ori A 


8.9 


0.5 


5.5 


0.8 


0.078 ± 0.005 


258.5 


83.9 


1672.6 


7860 ± 40 


TP 


A B 


3He 


G213.71-12.6 


MONR 2 


9.6 


1.3 


1.2 


0.5 


0.070 ± 0.002 


5.7 


40.5 


1567.9 


8570 ± 90 


SP 


A D 


3 He 


G220.508-2.8 


S291 


16.6 


9.2 


9.4 


25.2 


0.070 ± 0.002 


2.0 


0.2 


18.0 


21810 ± 4670 


SP 


C E 


3 He 


G223.7-1.9 




9.7 


1.6 


14.8 


6.7 


0.070 ± 0.002 


6.1 


0.3 


34.2 


9910 ± 1770 


SP 


C E 


3 He 


G227.79-0.12 


S298/NGC 2359 


12.4 


5.0 


10.6 


15.3 


0.099 ± 0.050 


3.5 


0.3 


25.4 


13560 ± 770 


SP 


A E 


•^He 


G231. 48-4.40 


RCW 6 


12.9 


5.8 


9.0 


15.1 


0.087 ± 0.020 


3.7 


0.4 


28.9 


9740 ± 260 


SP 


B D 


3 He 


G243. 16+0.37 


S311 


12.0 


5.4 


10.2 


16.0 


0.080 ± 0.009 


14.5 


1.4 


49.9 


10220 ±110 


SP 


B B 


3 He 


G345.03+1.54 




7.1 




10.8 




0.070 ± 0.002 


7.7 


0.6 




6850 ± 350 


SP 


D D 


Cll 


G345.23+1.03 




7.3 




4.4 




0.030 ± 0.002 


8.8 


4.4 




7590 ± 120 


SP 


B E 


Cii 


G345.31+1.47 




7.1 




7.3 




0.059 ± 0.004 


13.4 


2.4 




8530 ± 640 


SP 


B C 


Cii 


G345.40+1.41 




7.0 








0.059 ± 0.004 








... ± ... 




B D 


Cii 


G345.43-0.94 




6.2 


14.0 


3.2 


13.0 


0.059 ± 0.002 


27.5 


26.5 


228.4 


6960 ± 50 


SP 


A A 


Cii 


G345.54+0.10 




9.4 


17.5 


3.9 


19.9 


0.070 ± 0.002 


1.6 


1.0 


40.3 


12320 ± 970 


SP 


C E 


Cii 


G347.96-0.44 




2.5 




5.7 




0.070 ± 0.002 


2.3 


0.7 




5930 ± 200 


SP 


B D 


Cii 


G348.23-0.98 




6.3 


14.3 


3.3 


13.6 


0.182 ± 0.016 


8.3 


7.6 


118.6 


6610 ± 100 


SP 


A C 


Cii 


G348.72-1.03 




6.9 




4.3 




0.062 ± 0.004 


32.0 


17.3 




7150 ± 90 


SP 


B B 


Cii 


G350. 13+0.09 




2.6 


6.2 


3.1 


5.5 


0.072 ± 0.005 


6.1 


6.4 


171.7 


6710 ± 120 


SP 


B D 


Cii 


G351.063+0.6 




7.9 


0.6 


8.0 


1.4 


0.070 ± 0.002 


26.2 


4.0 


283.3 


10570 ± 340 


SP 


C D 


Cii 



Table 1 — Continued 



Source 


Name 


«gal 






U 


n\ xie )/n[^ti } 


c 


-iC 






iVioae 


i^r s 


Survey 






(kpc) 


(kpcj 


\ ) 










(cm~3) 










G351.192+0.7 




8.8 


17.1 


7.4 


36.9 


0.049 ± 0.003 


30.6 


5.4 


59.3 


5610 ± 20 


SP 


B A 


Cii 


G351. 201+0.4 




7.2 


1.4 


6.8 


2.7 


0.081 ± 0.007 


15.9 


3.4 


178.5 


6650 ± 70 


SP 


B C 


Cii 


0351.246+0.7'^ 


NGC 6334 A 


8.8 


17.1 


5.6 


27.7 


0.052 ± 0.005 


36.2 


11.4 


106.1 


8560 ± 70 


TP 


C B 


3 He 


G351.358+0.6 


S8 


7.9 


0.6 


7.0 


1.3 


0.067 ± 0.002 


54.4 


10.9 


460.7 


6840 ± 40 


SP 


C A 


Cii 


G351.368+0.7<1 


NGC 6334 D 


7.9 


16.2 


7.5 


35.2 


0.058 ± 0.004 


61.8 


10.9 


94.0 


9700 ± 90 


TP 


C A 


3He 


G351.467-0.4 




5.1 


3.4 


3.5 


3.4 


0.070 ± 0.002 


2.8 


2.3 


133.2 


7460 ± 120 


SP 


B C 


Cii 


G351.613-1.2 




6.0 


14.3 


2.9 


12.2 


0.088 ± 0.005 


23.2 


26.3 


238.2 


7620 ± 30 


SP 


A A 


Cii 


G351.64-1.26 




6.1 




3.1 




0.083 ± 0.004 


21.0 


21.5 




6490 ±210 


SP 


B B 


Cii 


G351.69-1.17 




6.4 


14.7 


5.1 


22.0 


0.077 ± 0.007 


11.8 


4.4 


72.2 


7560 ± 90 


SP 


B B 


Cii 


G352. 61-0.17 




2.0 


6.7 


4.1 


8.0 


0.070 ± 0.002 


2.6 


1.5 


70.3 


7560 ± 240 


SP 


B C 


Cii 


G353. 035+0.7 




6.8 


1.7 


3.7 


1.8 


0.116 ± 0.008 


7.1 


5.2 


260.2 


5630 ± 40 


SP 


B C 


Cii 


G353.08+0.36 




7.9 


0.6 


9.5 


1.8 


0.078 ± 0.009 


19.4 


2.1 


168.3 


5390 ± 100 


SP 


E C 


Cii 


G353. 14+0.66 




7.6 


16.0 


7.2 


33.5 


0.096 ± 0.006 


56.1 


10.6 


88.7 


6250 ± 30 


SP 


B B 


Cii 


G353.186+0.8 


Sll 


7.6 


0.9 


5.2 


1.3 


0.107 ± 0.004 


66.3 


23.7 


682.4 


7100 ± 40 


SP 


B B 


Cii 


G353.21+0.91 




7.7 




5.0 




0.099 ± 0.004 


62.1 


23.8 




6900 ± 50 


SP 


B B 


Cii 


G353.398-0.3 




5.2 


3.3 


2.6 


2.5 


0.050 ± 0.006 


9.3 


13.2 


376.7 


8480 ± 60 


SP 


A C 


Cii 


G353.43-0.37 




5.3 


3.3 


2.7 


2.5 


0.064 ± 0.011 


8.6 


11.8 


350.7 


7750 ±110 


SP 


B C 


Cii 



''The electron temperatures were rounded of to the nearest 10 K. 

''Mode refers to the continuum observing mode, where TP is for total power and SP is for switched power observations. 

"^QFs corresponds to the continuum and line quality factors, respectively. 

<^M16, M17, NGC 6334, Rosette and S209 were observed at different positions in the ^He survey. 
"^Observed in both C ii and ^He surveys. 



-31 - 



Table 2. Average Differences in Electron Temperature Determinations^ . 



Reference 


N 


100 < (us-ref)/us> 


100 < abs(us - ref)/us > 


Observation 


Method 


PTR78 


2 


3.6 


12.9 


optical 


[O in] A4363,5007 


WBW79 


6 


3.6 


4.9 


RRL 


(H 66 a)/Tc 


TMP80 


7 


2.9 


4.7 


RRL 


ri(H66a)/Tc 


TMP80 


5 


2.7 


11.0 


RRL 


ri(H 76 a)/Tc 


GR83 


3 


9.3 


11.1 


RRL 


Ti (H 125 a)/Tc 


WWB83 


30 


10.9 


14.4 


RRL 


Ti (H 76 a) /Tc 


SMNDP83 


17 


13.3 


14.2 


RRL 


Ti (H 109 a)/Tc 


ACA96 


4 


4.5 


15.8 


RRL 


Tl (H 66 a, 76 a, 93 a) /Tc 


BBRW99 


16 


-0.2 


6.7 


RRL 


Ti (H 91 a)/Tc 


DPCCOO 


2 


9.8 


9.8 


optical 


[Oiii] A4363,5007 


OCAR02 


2 


5.3 


8.7 


radio 


continuum 


EGPPRR05 


3 


5.1 


5.1 


optical 


[Oin],[S in], [Arm] 


EGPPRR05 


3 


-1.2 


5.5 


optical 


[On], [Nil] 



''QF = A, B and C, without the outUer G49.582-0.38. 



References. — ACA96 - AfHerbach et al. (1996); BBRW99 - Balser et al. (1999); DPCCOO - Deharveng et al. (2000); 
EGPPRR05 - Estcban ct al. (2005); GR83 - Garay & Rodriguez (1983); OCAR02 - Omar ct al. (2002); PTR78 - Pcimbcrt 
et al. (1978); SMNDP83 - Shaver et al. (1983); TMP80 - Thum et al. (1980); WBW79 - Wilson et al. (1979); WWB83 - 
Wink et al. (1983). 
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Table 3. T* Gradient Fits: T* = ai + oai^gai 



Sample 


ai(K) 


a2(Kkpc-l) 


r 




N 


Ai?gal(kpc) 


A=- 


5840 ± 410 


281 ± 49 


0.49 


1690 


109 


0.0-20.0 


Bb 


5660 ± 430 


286 ± 56 


0.50 


1430 


78 


0.0-17.0 


C<= 


5960 ± 380 


261 ± 50 


0.56 


1199 


64 


0.0-17.0 


Dd 


5480 ± 410 


312 ± 55 


0.56 


1053 


71 


2.0-13.0 


E"' 


5780 ± 350 


287 ± 46 


0.59 


1173 


76 


0.0-17.0 


pf 


5300 ± 350 


340 ± 45 


0.67 


1059 


73 


2.0-17.0 


-Rgal < Rq 


A* 


6110 ± 690 


222 ± 110 


0.23 


1675 


75 


0.0-8.5 


Bb 


6750 ± 660 


89 ± 106 


0.11 


1430 


62 


0.0-8.5 


0= 


6750 ± 550 


124 ± 90 


0.19 


1139 


52 


0.0-8.5 


Dd 


5730 ± 600 


268 ± 96 


0.35 


1007 


57 


2.0-8.5 




6540 ± 520 


153 ± 85 


0.23 


1130 


60 


0.0-8.5 


pf 


5730 ± 600 


268 ± 96 


0.35 


1007 


57 


2.0-8.5 


-Rgal > Rf) 




7190 ± 1370 


182 ± 115 


0.27 


1733 


34 


8.5-20.0 


Bb 


4590 ± 1510 


404 ± 130 


0.64 


1279 


16 


8.5-17.0 


c<= 


3780 ± 1910 


453 ± 157 


0.67 


1357 


12 


8.5-17.0 




5220 ± 2550 


342 ± 239 


0.38 


1305 


14 


8.5-13.0 


E« 


4590 ± 1510 


404 ± 130 


0.64 


1279 


16 


8.5-17.0 


pf 


4590 ± 1510 


404 ± 130 


0.64 


1279 


16 


8.5-17.0 



''AH QPs: 106 sources plus 6 additional pointings within 5 of those sources. 
S291 and S132 were not included due to their large errors; we have no continuum 
for G345.40+1.41. 

bQF = A, B and C. 

'^Sample B excluding sources calculated with only TP continuum, 
d Sample B excluding outlying values of Te* flagged in Figure 3. 
"Sample B excluding the outliers G49.582-0.38 and G5.956-1.265. 
^Sample E, but excluding Sgr B2, Gl.13-0.1 and G5.899-0.427 
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Table 4. T* Gradient Fits: T* = ai + a2-Rgai for Different Azimutfi Intervals^ . 



Azimuth Range 


ai(K) 


a2(Kkpc-i) 


r 




N 


Aii (kpc) 


330°-30° 
















5730 ± 480 


294 ± 53 


0.73 


997 


29 


2.0-16.9 




6070 ± 1070 


233 ± 157 


0.33 


970 


20 


3.0-9.0 


30° -90° 
















4540 ± 650 


493 ± 95 


0.79 


1010 


18 


3.7-12.3 




4830 ± 950 


432 ± 168 


0.58 


990 


15 


3.7-9.0 


90° -150° 
















4320 ± 2050 


507 ± 324 


0.48 


1320 


10 


4.4-8.5 




4320 ± 2050 


507 ± 324 


0.48 


1320 


10 


4.4-8.5 


150°-215° 
















8270 ± 1360 


-92 ± 200 


-0.13 


1481 


14 


0.4-8.8 




8500 ± 2910 


-120 ± 410 


-0.09 


1546 


13 


3.0-8.8 



''QF = A, B and C, without the outliers G49.582-0.38 and G5.956-1.265. 
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Table 5. Comparison of Electron Temperature Gradient Determinations 



Reference 


Ail 


dTe/dR 


N 


Observation 


Method 


PTR78 


8.4-13.9 (10.0) 


11000 


5 


optical 


[Oiii] A4363,5007 


LRC79 


5.0-13.0 (10.0) 


390±70 


20 


RRL 


Ti (H 86 a)/Tc 


LRC79 


5.0-13.0 (10.0) 


250±150 


20 


RRL 


Ti (H 108 /3)/Tc 


WPZ79 


5.0-12.5 (10.0) 


400±100 


18 


RRL 


Ti,(H65a)/Tc 


DWBW80 


4.0-09.0 (10.0) 


340±90 


US'* 


RRL 


Ti (H 110 a)/Tc 


MN81 


5.0-11.0 (10.0) 


420 


23 


RRL 


Tl{R 7&a)/Tc 


GR83 


3.6-12.6 (10.0) 


440±50 


23 


RRL 


ri,(H125«)/Tc 


SMNDP83 


3.5-13.7 (10.0) 


433±40 


67 


RRL 


r£,(H 109 a, 76 a)/Tc 


WWB83 


4.0-17.0 (10.0) 


270 


84 


RRL 


Tl (H 76 a)/Tc 


ACC85 


4.0-12.0 (10.0) 


310 


27 


RRL 


Ti, (H 166 «)/Tc 


ACA96 


4.0-11.0 (8.5) 


320ib64 


17 


UC RRL 


T^in 66 «, 76 Q, 93 «)/Tc. 


DPCCOO 


6.6-14.8 (8.5) 


372±38 


6" 


optical, RRL 


[Oiii] A4363,5007, Tl/Tc 


This work 


0.0-17.0 (8.5) 


287 ± 46 


76 


RRL 


r£,(H 91 a)/Tc 


This work (-GC) 


2.0-17.0 (8.5) 


340 ± 45 


73 


RRL 


ri(H91 a)/Tc 



^Objects averaged in 6 groups corresponding to increments of 1 kpc in Galactocentric radius. 

^^£([0111]) in 6 H II regions plus mean temperatures obtained, in 1 kpc bins, from the following 
RRLs studies: Mezger et al. (1979), Shaver et al. (1983), Wink et al. (1983) and Caswell & Haynes 
(1987). 



References. — ACA96 - Afflerbach et al. 1996; ACC85 - Azcarate et al. 1985; DPCCOO - Deharveng 
et al. 2000; DWBW80 - Downes et al. 1980; GR83 - Garay & Rodriguez 1983; LRC79 - Lichten et al. 
1979; MN81 - McGee & Newton 1981; PTR78 - Peimbert et al. 1978; SMNDP83 - Shaver et al. 1983; 
WPZ79 - Wilson et al. 1979; WWB83 - Wink et al. 1983. 
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LTE Electron Temperaturejiterature (K) 



Fig. 1. — Nebular LTE electron temperatures derived here compared with those from Shaver 
et al. (1983) (circles) and Wink et al. (1983) (triangles) for objects in common. Filled 
symbols denote our highest quality data (QFs A, B, Sz C); open symbols flag our poorer 
quahty data. The locus of Te(us) = T'e(hterature) is shown by the solid hne. 




Fig. 2. — Distribution of the H ll region sample projected onto the Galactic plane plotted as a 
function of Galactocentric radius and azimuth. Different symbols identify the ^He (triangles) 
and C II survey (circles) nebulae. Filled/open symbols have the same meaning as in Figure 1. 
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Fig. 3. — Nebular LTE electron temperature plotted as a function of the Galactocentric 
distance for the best data (QFs C or better — Sample B of Table 3). Least squares linear 
fits to the gradient for the entire sample and, separately, for nebulae located inside and 
outside the solar orbit are shown. The gradient is flatter for H ll regions in the inner Milky 
Way. 




-38- 



14000 




4 6 8 10 12 14 16 18 20 
Galactocentric Distance (kpc) 



Fig. 4. — Repeat of Figure 3 showing the quahty factor of the continuum data. Filled 
symbols denote T* temperatures derived from switched power continuum measurements; 
open symbol sources used total power measurements. 
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14000 I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




4000 I I ' I I I ' I I I ' I I I ' I I I ' I I I 

2 4 6 8 10 12 14 16 18 20 

Galactocentric Distance (kpc) 

Fig. 5. — "Smoothed" electron temperatures plotted as a function of Galactocentric distance, 
RgaX- Symbols with error bars are a 10 point running mean of T* and i?gai drawn at intervals 
of 4 points along -Rgai- Horizontal line segments show the mean T* in each 1 kpc wide i?gai 
interval. These are offset by +2,000 K for clarity. Both smoothing algorithms suggest a 
slightly smaller gradient in the inner Galactic disk. 
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-60 -40 -20 20 40 60 80 
100 [Te(LTE) - Te(sample E)] / Te(LTE) 

Fig. 6. — Histogram of the fractional deviation of the nebular electron temperature from the 
best temperature gradient model fit to sample E. The deviations are well fitted by a Gaus- 
sian distribution whose dispersion (sigma) corresponds to an intrinsic electron temperature 
fluctuation of ~ 1,100 K at any Galactic radius. 
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Fig. 7. — Nebular electron temperature plotted as a function of Galactocentric radius for 
four distinct ranges of Galactocentric azimuth. The analysis is compromised because our 
Hll region sample is not uniformly distributed in the Milky Way's disk (see Figure 2). 
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Fig. 8. — Comparison of our Hll region electron temperatures with those from Shaver 
et ah (1983). The Shaver et ah nebulae were observed with the H 76 a (triangles) and 
H 109 a (circles) recombination lines. Their sample has a large number of objects in the 
interval i?gai = 2-6 kpc. Their H 109 a T* values in this zone are systematically lower than 
ours which may contribute to the very steep electron temperature gradient they derive. 
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Fig. 9. — The Galactic 0/H abundance gradient derived from our nebular electron temper- 
atures using the relation between 0/H and Tg of Shaver et al. (1983). Symbols have the 
same meaning as in Figure 3. The sohd hne is the Eqn. 6 least squares fit. 



